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Abstract 


We have treated in detail the theory of gamma-ray line production in 
solar flares. The strongest line, both predicted theoretically and detected 
observationally at 2.2 MeV, is due to neutron capture by protons in the 
photosphere. The neutrons are produced in nuclear reactions of flare 
accelerated particles which also produce positrons and prompt nuclear 
gamma rays. From the comparison of the observed and calculated intensities 
of the lines at 4.4 or 6.1 MeV to that of the 2,2 MeV line jt is possible 
to deduce the spectrura of accelerated nuclei in the flare region; and 
from the absolute intensities of these lines it is possible to obtain 
the total number of accelerated nuclei at the Sun. The study of the 2.2 MeV 
line also gives information on the amount of He 3 in the photosphere. The 
study of the line at 0.51 MeV resulting from positron annihilation complements 
the data obtained from the other lines; in addition it gives information 
on the temperature and density in the annihilation region and on the 
anisotropy of the accelerated electron beam which produces continuum gamma 
rays at energies greater than about 1 MeV. 
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1. INTRODUCTION 

Measurements of accel erated charged particles near the Earth clearly 
indicate that such particles are produced in great profusion in solar 
flares. These particles consist of both electrons and nuclei; but until 
the advent of solar gamma-ray astronomy, observations in the radio and 
x-ray bands had revealed only the existence of the electronic component 
in the flare region itself. 

In a previous paper (Lingenf elter and Ramaty 1957) we treated in 
considerable detail the nuclear reactions produced by accelerated charged 
particles in solar flares and we showed that in large flares these 
reactions produce detectable lines in the gamma-ray region. We found 
that the strongest lines should be at 0.5, 2.2, 4.4, and 6.1 MeV resulting 
from position annihilation, neutron capture on, hydrogen, and deexcitation 
of excited states in C 12 and 0 16 , respectively. 

The recent observations by Chupp et al. (1973) of the first gamma- 
ray lines from solar flares confirm these predictions. During the flash 
phase of the 1972, August 4 flare all of these lines were observed with 
relative intensities essentially consistent with our calculations. 

Since these observations became available, several additional studies 
on gamma-ray line production in solar flares have been undertaken. 

Ramaty and Lingenf elter (1973a) have investigated the consistency of 
the observations with the theory of nuclear reactions in flares; these 
authors also considered the effects of positronium formation and neutron 
propagation in the solar atmosphere (Ramaty and Lingenfelter 1973b); 
Reppin et al. (1973) treated the time dependence of the 2.2. MeV line; 
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Wang and Ramaty (1974) have done a detailed calculation on neutron propa- 
gation and 2.2 MeV line formation and they pointed out the importance 
of photospheric He 3 as a nonradiative sink for the neutrons; Kozlovsky 
and Ramaty (1974a) have pointed out that aa reactions produce the Li 7 and 
Be 7 lines at 478 keV and 431 keV which could be observable from flares; and 
Ramaty and Kozlovsky (1974) evaluated in detail the production of H 2 , H 3 
and He 3 in flares and they attempted to deduce the number of protons released 
from the flare of 1972, August 4 by combining the He 3 and gamma-ray obser- 
vations. 

In the present paper we wish to summarize the above material and to 
present updated calculations on the production of gamma-ray lines in 
solar flares. In Section 2 we define the interaction models that we 
use in our calculations; in Section 3 we consider neutron production and 
2.2 MeV line formation; in Section 4 we consider the production of prompt 
gamma- ray lines with special emphasis on the 4.4 MeV and 6.1 MeV lines for 
which observational data exists; in Section 5 we compare the results of 
Sections 3 and 4 with data for the 1972, August 4 flare and we deduce 
the number and spectrum of accelerated particles at the Sun; in Section 6 
we treat problems concerning the formation of the 0.51 MeV line, and we 


summarize our results in Section 7. 
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2. INTERACTION MODELS 

We consider two limiting interaction models (e.g. Ramaty and 
Lingenfelter 1973a)' A thin-target model in which the spectrum of 
accelerated particles is not modified during the time in which the 
nuclear interactions take place, and a thick-fa-rget model in which 
the accelerated particles move from the flare region downward into the 
Sun, undergoing nuclear interactions as they slow down in the solar atmos- 
phere. In the thin- target model it is assumed that either the total 
path length traversed by the particles at the Sun is small in comparison 
with their interaction length, or that the particle energy Joss from 
ionization and nuclear interactions is just balanced by energy gains 
from acceleration. 

For the composition of the ambient solar atmosphere we use the abundances 
given by Cameron (1973). For the accelerated particle populations in both 
the thin and thick-target models we consider power-law and exponen- 
tial spectra. In the thin-target model these are 

N ± (E) = k^ -5 , (1) 

and 

/ 

N t (P) - K ± exp(-P/P 0 ), (2) 

respectively. Here Nji(E) and N^(P) are the instantaneous numbers of 
accelerated pai tides of kind i in the interaction region per unit 
energy per nucleon, E, or unit rigidity, P; and k^ are constants 
determined by normalizing the N^ ' s to 1 proton of energy greater than 
30 MeV and by using the composition of the ambient solar atmosphere; 
and s and P Q are, respectively, the spectra 1 - index and characteristic 
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rigidity assumed to be the same for all accelerated particle components. 
In the thick- target model we use expression similar to equations (1) 
and (2), but we replace the instantaneous numbers by total numbers, 

Nj, such that N^(E) and N.. (P) are the total number of accelerated charged 
F articles per unit energy per nucleon or unit rigidity that are released 
from the flare region downward into the Sun. As with the instantaneous 
fluxes in the thin-target model, the normalizations of the ISL 's are 
determined by using the composition of the ambient solar atmosphere and 
1 proton of energy greater than 30 MeV. 

In the thin- target model, the production rate of secondaries from a 
particular reaction is given by 

q i = n i dE v E) ^ ct i ( e) (3) 

where n-^ is the number density of target atoms in the solar atmosphere, 
cjB is particle velocity, and aj[(E) is the cross section as a function of 
energy per nucleon. The units of are secondary particles per second. 

In the thick-target model, the total production of secondary particles 
for a given reaction is 

Q i * Tii f dp < e '> dE % °i< E > • < 4) 

O O 

where J| is the number of target nuclei per gram of solar material, and 
dE/dx is the stopping power of the primary particles in solar material 
due to both Coulomb and nuclear collisions. The quantities are the 
total number of secondaries . Expressions similar to equations (3) and 
(4) can be written down for rigidity spectra. 
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By inverting the order of integration, equation (4) can be written 


as 


^ dx 

Q x = \ Jo dE di 


a,(E) f dE' N. 0 ? )• 
1 E 1 


(5) 


Because the inner integral is just the integral spectrum of N^E), 
equation (5) can be further simplified, 

«i -HiJT dE i°i <E) 5 i <>«• 


( 6 ) 


In our subsequent treatment of secondary particle production, we shall 
use equations (3) and (6) or their equivalents for rigidity spectra. 



-7- 


3. NEUTRON AND 2.23 MeV GAMMA RAY PRODUCTION 

Neutron production by accelerated charged particles was treated by 
Lingenfelter et al. (1963) and Lingenfelter and Ramaty (1967)- Recently, 
Kozlovsky (private communication) has updated the cross sections used by 
these authors and has added additional data at low energies including data on 
alpha particle induced reactions. The neutron production cross sections 
are shown in Figure 1. Here pp, pa*, aar, pCNO and aCNO indicate neutron 
production in proton-hydrogen, proton-helium, alpha particle-helium, proton- 
heavy nuclei, and a particle-heavy nuclei reactions, respectively. The 
latter two cross sections are the neutron production cross sections for 
all nuclei with A^12, normalized to one such nucleus by using the 
el^mertal and isotopic abundances of Cameron (1973). 

The instantaneous neutron production rates in the thin-target model 

for power- law and exponential spectra are shown in Figures 2 and 3, 

respectively. The various production modes are: pp (proton-hydrogen), 

pa (proton-helium), ap (a particle-hydrogen), aa (a particle-helium), 

pCNO (proton-heavy nculei) , CNOp (heavy nuclei-hydrogen) , aCNC ( a particle- 

heavy nuclei) , and CNQa (heavy nuclei-helium) . As can be seen for flatter 

spectra (smaller values of s or larger values of P^) the neutrons are 

produced mainly in pa, pp and ap reactions. For power-law spectra (Figure 2) 

neutron production at large values of s is mainly due to aCNO and CNQa 

reactions; the contribution of pCNO and CNOp reactions is small at all 

values of T; and ota reactions make a major contribution around s*4. For 

exponential spectra, (Figure 3) almost all of the neutrons are produced in 

pa reactions at most values of P . In this case the relative contributions 

o 

of reactions induced by a-particles and heavy nuclei are lower than for 



power-law spectra. Because particles with Z^2 have larger rigid:* ti.es 
than protons, it follows from equation (2) ^hat their fluxes relative to 
the proton flux at the same energy per nucleon is lower than in the power- 
law case given by equation (1) . Since the nuclear cross sections are the 
same for the direct and inverse reactions at the same energy per nucleon, 
for particle spectra which are exponential in rigidity more neutrons are 
produced by a proton induced reaction than by the corresponding inverse process. 

The total neutron production rates in the thin-target model, q, and tne 
total neutron yields in the thick-target model, Q, are shown in Figure 4 
for power-law and exponential spectra. In the thin- target model, for flat 
primary spectra the neutron production rates are about the same for the 
exponential and power -law cases. This result is due simply to the fact that 
for such spectra most of the neutrons are produced in pa reactions with 
effective threshold around 30 MeV/nucleon and all of the assumed spectra 
are normalized of 1 proton above this energy. But because a power law 
spectrum contains a much larger number of low-energy particles than an 
exponential spectrum with this same normalization, the steep power law 
spectrum can yield orders of magnitude more neutrons from orCNO and CNQa 
reactions which have thresholds more than an order of magnitude below 30 
MeV/nucleon. 

Similar effects are evident also in the thick- target model except that 
Here the relative contribution of the low-energy particles in general is 
diminished because of their shorter ran 0 ^. 

Having considei*ea !~he production of neutrons, let us now discuss their 
propagation and the e^ suing gamma-ray line production. Wang and Ramaty 
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(1974) considered in detail the effects of neutron propagation in the 
solar atmosphere on the production of gamma rays by the reaction 

n + p d + y . (7) 

In their treatment a distribution of neutrons was released in the 
chromosphere or corona, and the path of each neutron after its release 
was followed by a computer Monte-Carlo simulation. If the neutrons are 
released above the photosphere, any initially upward moving neutron 
escapes from the Sun. Some of the downward moving neutrons can also 
escape after being backscattered elastically by ambient protons, but most 
of these neutrons either are captured or decay at the Sun. Because the 
probability for elastic scattering is much larger than the capture 
probability, the majority of the neutrons are thermaJ ized before they 
get captured. Since the thermal speed in the photosphere (where most 
of the captures take place) is much smaller than the speed of light, 
the gamma-rays from reaction (7) are essentially all at 2.2 Me u and the 
Doppler-broadened width of this line is negligible. 

The bulk of neutrons at the Sun are captured either on H or on He 3 . 
Whereas capture on H yields a 2.2 MeV photon, capture on He 3 proceeds via 
the radiationless transition 

n + He 3 -4 H 3 + p, (8) 

and hence produces no photons. The cross sections for reactions (7) and 
(8) are 2.2 x 1CT 30 g"* 1 cm 2 and 3.7 x lCT^g^cnF, respectively, where g is 
the velocity of the neutron (for details see Wang and Ra. aty 1974) . Thus 
if the He 3 /H ratio '»n the photosphere is ~ 5 x 10 5 comparable to that 
observed in the solar wind, nearly equal numbers of neutrons are captured 
on He 3 as on H. 



- 10 - 


The results of the Monte-Carlo calculations of Wang and Ramaty (1974) 
are presented in Figures 5 and 6 for two assumptions on the phozospheric 
He 3 abundance: He 3 /II = 0 and He"' /R = 5 x 10” 5 . In these calculations 

an isotropic distribution of monoenergetic neutrons of energy E a is 
released above the photosphere. The solid lines are the probabilities 
for the various indicated processes. As can be seen, the capture and 
loss probabilities increase with increasing energy, because higher energy 
neutrons penetrate deeper into the photosphere. This reduces their escape 
probability and leads to a shorter capture time, thereby reducing the 
decay probability Whan He 3 /H - 5 x 10~ 5 , the probability for loss on 
He 3 almost equals the capture probability on protons. The escape probability 
is greater than 0.5, because all initially upward moving neutrons escape 
from the Sun. Noce that the sum of all probabilities equals 1. 

The dashed lines in Figures 5 and 6 are photon yields per neutron, 
f(0,E n ), for various neutron energies, E^, and angles, 0, between the 
earth-sun line and the vertical to the solar surface. The function f is 
defined such that for an average neutron production rate, q, the average 
2.2 MeV photon flix at Earth is 

* (2.2 MeV) = qf/(4rrF 3 ) (9) 

where R — 1 A.U. 

At low neutron energies and e near zero, f is close to the capture 
probability on protons. This means that gamna rays from low-energy 
neutrons observed close to the vertical escape essentially unattenuated 
from the Sun, At higher energies and at larger angles, however, there 
is significant attenuation of the gamma rays due to Compton scattering 
in the photosphere. Even though f does depend on E^, for flares 
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sufficiently close to longitude and latitude zero on the Sun and 
neutron energies between about 1 and 100 MeV, we can approximate it 
by a constant. Most of the neutrons have energies in this - ange 
(Lingenfelter and Ramaty 1967). Thus, for He 3 /H=~ 5 x 10“ 5 we use 
f =- 0.12, and for He 3 /H^ 0, we take f ^ 0.2. Note that these approxi- 
mations are quite valid for the flare of 1972, August 4, since its solar 
longitude and latitude where E08 and N14. 

It should be noted that equation (9) is valid for the average neutron 
flux only, because the instantaneous 2.2 MeV flux lags behind the instan- 
taneous neutron production rate. This lag is almost entirely due to 
the finite neutron capture time in the photosphere. Wang and Ramaty (1974) 
have investigated this effect, and some of their results are given in 
Table 1. Here <n> is the most probable density in the photosphere where 
the captures take place, t is the mean capture time, and T , is the neutron 

^ Cl 

decay mean life. In terms of the parameter A. = r c -1 + r > the time 
profile of the 2.2 MeV photon flux from a monoenergetic burst of neutrons 
released at t Q can be approximated (Wang and Ramaty 1974) by 


(2.2 MeV) * expf-X (t-t Q ) ] . 


( 10 ) 



- 12 - 


4 . PROMPT GAMMA-RAY LINES PRODUCTION 

The various prompt gamma-ray lines that can be produced in solar 
flares together with their production mechanisms are listed in Table 2. 

These lines have already been discussed by Ramaty and Lingenfelter 
(1973a) except for the Li 7 and Be 7 lines at 478 keV and 431 keV. The 
possibility of producing these lines in ckx reactions has been pointed out 
recently by Kozlovsky and Ramaty (1974a) 

The cross sections for the reactions C l2 (p,p') C 12 * 4> 4si , C 12 (osc/ ) C 12 ***• 43 , 

0 16 (p,p / ) O 16 * 6 * 14 , and O lG (a,c/) 0 16 * 6 * 14 are given in Figure 7. 

The cross section for the proton induced reactions were summarized by 
Lingenfelter and Ramaty (1967), and the cross sections for the cr-particle 
induced reactions are from Kozlovsky (private communication). In addition, 
we also consider the reaction p+O l6 -*C 12 * 4# 4o +. . • which was discussed 
previously (Ramaty and Lingenfelter 1973a). By using these cross sections 
and the interaction models discussed above, we can calculate the production 
rates of C 12 * and 0 16 * and the resultant prompt photons at 4.43 MeV and 
6.14 MeV. We must distinguish, however, between reactions induced by 
accelerated protons or o-particles, and reactions induced by accelerated 
heavy nuclei. For the former, the Doppler widths of the lines are small 
in comparison with available instrumental resolutions (about 130 keV) . 

But for the latter the lines are significantly broadened by the motion 
of the excited fast nucleus which has lost little kinetic energy in die 
interaction. Because these lines are so broad that they cannot be resolved 
from the background with presently available instrumentation, in our 
treatment we consider the intensities of the 4.43 MeV and 6.14 MeV lines 
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from proton and ^-particle induced reactions only. 

The production rates of C 12 * and 0 16 * in the thin-target model 
with power-law spectra are shown in Figure 8. We see that for flat spectra, 
the lines at 4,43 and 6,14 MeV are produced mainly by proton induced reactions, 
whereas for steep spectra (large values of s) , the contributions of the 
a-par tides becomes important. For exponential spectra, pro ton- induced 
reactions are the principle source of the excited states at all values of 

v 

The cross section for Li 7 production in oo' reactions is shewn in 
Figure S (Kozlovsky and Ramaty 1974b), The c^oss section for Li 7 * produc- 
tion is about half of the total Li 7 production independent of energy. 

Similarly, the cross section for Be 7 * should also be about half of the 
total Be 7 production; and because Be 7 and Li 7 are produced by mirror 
reactions their production cross sections should be about equal, even though 
no data for Be 7 production in aor reactions is avialable. 

Using these cross sections we find that the intensities of the lines 
at 478 keV and 431 keV are approximately the same as the intensity of 
the 4.43 MeV line. Their Doppler width, however, are about 30 keV 
(Kozlovsky and Ramaty 1974a); therefore, they should not be observed 
individually but rather as a broad spectral feature. 

The cross section for the other lines of Table 2 and their intensities 
were discussed by Ramaty and Lingenfelter (1973a). The intensities of 
these lines are significantly lower than the intensity of the 4.43 MeV 
line for all interaction models and spectra . 
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5. ACCELERATED PARTICLES AT THE SUN AND IN THE INTERPLANETARY MEDIUM. 

The energy spectra of the accelerated particles at the Sun can be 
deduced by comparing the calculated and observed ratios of the intensities 
of the strongest prompt line at 4.43 MeV to that of the neutron capture 
line at 2.23 MeV. This ratio for both the thin and thick- target models 
is shown in Figure 10 as a function of s, for power-law spectra, and P c , 
for exponential spectra* These were obtained from the calculations of 
Sections 3 and 4 using a photon yield, f, of 0.2. This yield corresponds 
to a photospheric ratio He 3 /H = 0. For He 3 /H * 5xlG ~ 5 the curves in 
Figure 10 should be raised by about a factor of 2. The 4.43 MeV line is 
expected to be somewhat stronger than the 6.14 MeV line. Since the same 
observational limits were reported for these two lines, t we normalize 
our calculations to the 4.43 MeV line, our results will also be consistent 
with the 6.14 MeV line within the uncertainties of the measurements. 

As can be seen from Figure 10 for exponential spectra in both the thin 
and thick- target models $4 43^2 23 decreases with increasing P Q . This 
results from the increase of the neutron production cross section with 
increasing energy as opposed to the decrease of the excitation cross 
sections (compare Figures 1 and 7). The same behavior can be seen for 
power-law spectra for s smaller than about 4.5 in the thin-target model, 
and s less than about 6 in the thick- target model. For larger values of s, 
$4 43^2 23 decreases with increasing s because these neutrons are 
produced mainly in aCNO and CNOa reactions which have lower thresholds than 
those for prompt gamma-ray production. 
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Let us compare now the calculations with the data. The observed 
(Chupp 1974) 43/^2 23 rat -*-° 0.11 +0.04 for the 1972, August 4 flare 

is also shown in Figure 10. As can be seen from this figure, assuming power- 
law spectra for the particles in the flare region, their spectral index 
s should lie between the values of 2 4^.2 deduced for the thin- target 
model and 3+0.3 for the thick- target model if there is no He 3 in the 
photosphere; or between 1.8 j0.2and 2.7 +0.2 for these models if the 
photospheric He 3 /H ratio is 5xl0“ 5 . Similarly, assuming exponential 
spectra, the implied Pq’s should lie between 110 +30MV for the thick- target 
model and 160 +35MV for the thin target model if He 3 /H = 0; or between 
180 +50MV and 230 +50MV for these models if He 3 /H = 5xl0~ 5 . 

Compart on of these implied particle spectra in the solar flare region 
with the proton spectrum observed in the interplanetary medium from the 
1972, August 4 flare is complicated by the possibility that the latter 
spectrum may have been significantly modified by acceleration in the 
interplanetary medium. Nonetheless, the proton spectrum obtained by 
Bertsch et al. (1974) from a rocket flight at 1916 UT on August 4, 1972 had 
a spectral index s=*2 in the 10 to 100 MeV region, in good agreement with our 
conclusions . 

Knowing the spectral index for the various models, we can now 
deduce t Ye total number of protons at the Sun. In the thick- target model, 
he ;ime integrated photon flux from the flare determines the total number 
of accelerated particles that interact and stop at the Sun. According 
to Chupp (1974), the 2.2 MeV intensity for the 1972, August 4 flare was 
about 0.3 photons cm“ 2 s* 1 . This flux is the average over the time 
interval of observation (0623 to 0633 UT) which gives a total flux of 
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about 180 photons cm -2 . Since the detector on 0S0-7 was eclipsed 

by the Earth before the termination of the gamma-ray event, this is a lower 

limit to the total flux from the flare. If we assume a duration of ^10 3 

seconds for the event, as indicated by the hard x-ray data (van Beek et al. 

1973), then the total flux of 2.2 MeV photons was ~300 cm*" 2 . From the 

neutron yield Q shown in Figure 4, we find that if f = 0.2, N^(>30MeV) - 

7X10 32 for s=3, and N (>30MeV) = 1.0 x 10 33 for P = 110 MV. The corres- 

P o 

ponding values for He 3 /H = 5xl0~ 5 are about the same, because the effect 
of a lower f is approximately cancelled by a higher neutron yield . Therefore, 
the total number of protons above 30 MeV released downward into the Sun 
in the thick-target model is about 10 33 independent of the spectral form 
or the He 3 /H ratio. 

For the thin- target model, the instantaneous gamma-ray observations 
determine the product of the ambient proton or hydrogen density, n^, and 
the instantaneous number of the accelerated particles in the interaction 
region. From the neutron yield q shown in Figure 4, the flux given above 
implies that if f - 0.2, n^Np (>30MeV) 5xl0 43 cnT 3 for s=2, and n^N^ 

(>30MeV) =- 10 44 cnf 3 for P ■ 160MV. As before, these values are not 

o 

effected by the photospheric He 3 /H ratio. 

To obtain an estimate of the number of protons released into the 
interplanetary medium, the information obtained from the gamma rays in 
the thin-target model can be combined with data on the path length traversed 
by the nuclei before their escape from the interaction region. Such 
information can be obtained from studies of deuterons and helium-3 nuclei 
from flares. The H 2 and He 3 observations from the 1972, August events 
(Webber et al. 1974), when compared with calculations of the production of 
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these isotopes in nuclear reactions (Ramaty and Kozlovsky 1974) imply that 
the amount of matter traversed by relativistic particles is about 1.5 g cm -a . 
This means that the product is about 3X10 13 cnT 3 s. where t : is the 

interaction time of the particles at the Sun. If tj is also interpreted 
as the escape time of the particles from the interaction region, then 
the protons were released into the interplanetary medium at an average 
rate N^(>30MeV) / 1^ varying from about 1.6 to 3X10 30 protons s -1 . The 
total number of protons released is the product of this rate and the 
acceleration time T. As indicated by the x-ray data (van Beek et al. 1973), 

T is about 10 3 seconds; therefore the total number of protons released is 
between about 1.6 to 3X10 33 protons. These numbers are larger by only 
about a factor of 2 to 3 than the number of protons released downward into 
the Sun in the thick-target model. 

The number of protons released from the flare should be compared with 
estimates of the total number of protons in the interplanetary medium as 
obtained from charged particle observations. Ramaty and Lingenfelter (1973a) 
have estimated this quantity by taking the observed density of protons 
near Earth and by multiplying it with a storage volume in the interplanetary 
medium which they took to be ~ lu 39 cm 3 . By using this volume and the peak 
proton flux greater than 30 MeV measured by Kohl et al. (1973), we get 
about 10 34 protons. This number is larger by a factor of 3 to 6 than our 
estimate for the number of protons released in the thin- target model. 

But since the measured peak proton flux could consist to a large degree of 
particles accelerated by shocks in the interplanetary medium, it appears 
that there is no real discrepancy between the number of protons released 
for the Sun as deduced from the gamma rays and the number observed in 
interplanetary space. 
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6. THE NATURE OF THE POSITRON ANNIHILATION RADIATION 

Positrons can be produced in solar flares from the decay of rr 
mesons and radioactive nuclei which result from nuclear reactions of 
accelerated protons and nuclei with the ambient medium. The possibility 
that a significant fraction of solar positrons can also be produced by 
accelerated electrons in e + -e pairs is being considered by Bai and 
Ramaty (1975). Annihilation of positrons can produce a gamma-ray line 
at 0.51 MeV. This line was observed by Chupp et al. (197 d) for both 
the 1972, August 4 and August 7 flares. 

The cross sections for the production of n + mesons and g + -emitting 

radioactive nuclei were given by Lingenfelter and Ramaty (1967). Kozlovsky 

(private communication) has recently updated the cross sections of the 3 - 

emitters. Using these cross sections we calculate from equation (3) whe 

instantaneous production rates of positron-emitters in the thin-target 

model for power law spectra. The results are shown in Figure 11. As can be 

seen, for flatter spectra the principal positron source is tt + mesons, but for 

steeper spectra, the pion contribution is negligible. This effect can also 

be seen in Figure 12, where the ratio of the pion yield to the total 

positron-emitter yield (including both pions and radioactive nuclei) is 

plotted as a function of s or P q for both interaction models. We see 

that, as for neutron production, the efficiency of tt + production is larger 

in the thick- target model than in the thin- target model. The ratios 
yield 

of the total/of positron-emitters to the neutron yield are shown in 
Figure 13 for the various interaction models. 

The intensity of the 0.51 MeV line, however, depends not only on 
the number of positron emitters produced, but also on the decay rate of the 
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po^itron emitters and on the annihilation rate of the positrons. The 
mean lives against radioactive decay of the various positron emitters 
are shown in Table 3. 

The annihilation rate of the positrons depends on the density, temperature 
and state of ionization of the ambient medium. As discussed previously 
(Stecker 1969, Ramaty and Lingenf elter , 1973b) in a low density and 
neutral medium, positrons nearly always annihilate via positronium forma- 
tion; only if the ambient density exceeds about 10 15 cm"* 3 do collisions destroy 
the 3 S state of positronium at such a rate that free annihilation becomes 
more important than positronium annihilation (Leventhal 1973). Positron 
annihilation from a bound state of positronium results in an asymmetric 
0.51 MeV line, since 75% of the time positronium annihilates from the 3 S 
state into 3 photons of energies less than 0.51 MeV, instead of 2 photons 
at precisely this energy. The 0.51 MeV line from the Sun, however, appears 
to be symmetric, and the observational upper limit on the fraction of 
3-photon annihilations is about 20% of the total (D . Forrest, private 
communication 1973). 

As we noted, three-photon annihilation could be reduced by collisions, 
but the required density of ^ 10 15 cm"* 3 is quite large. However, if the 
ambient medium is ionized, the rate of positronium formation is greatly 
reduced. In Figure 14 we show the rate of positronium formation and free 
annihilation per positron in a hydrogen plasma of unit density as a function 
of its temperature (C. Werntz and C. Crannell, private communication 1973). 

As can be seen, if the temperature is greater than about 10P°K, most 
of the positrons annihilate without forming positronium. This is not an 
unreasonable temperature for the annihilation region; however, from the 
observed upper limit on the width of the 0.51 MeV line, the temperature 
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in the annihilation region should be less than ~ 10 7 °K (Chupp 1974). 

According to the observations of the 1972, August 4 flare, the rise 
time of the 0.51 MeV line was within instrumental errors, similar to or 
perhaps even shorter than the rise time of the 2.23 MeV line (Chupp 1974). 

The latter was about 100 seconds, consistent with the expected lag between 
the production of the neutrons and the formation of the 2.2 MeV line 
(see Table 1) . Assuming that the time dependence of the production rate 
of the positron emitters is the same as that of the neutrons, the lag 
between the production of the positron emitters and the formation of the 
0.51 MeV line should not be longer than about 100 seconds. From Figure 14, 
this result implies that the density of the ambient medium in the annihil- 
ation region is at least 10 1S cmT 3 . Furthermore if the positrons are of 
nuclear origin, then from Table 3 if follows that they should mainly 
result from rr + mesons and short lived radioactive nuclei (0 1S , 0 14 and N 12 ). 

The observed average flux in the 0-51 MeV line for the flare of 1972, 
August 4 was about 0.06 photons cnT 2 s -1 (Chupp 1974). Hence the 

observed t>Q 52/^2 23 rat ^° was a ^ out 0*2. 

Let f 7 be the 0.51 MeV photon yield per positron defined in the 
same way as the 2.23 MeV photon yield per neutron, f, in equation (9). 

As discussed in Section 3, f ranges from about 0.1 to 0.2. The maximum 
value of f 7 is 2. But because part of the positrons can escape from the 
Sun before they annihilate, and, furthermore, a fraction of the positrons 
can be trapped at the Sun in low density regions where the annihilation 
time is long, f 7 can be considerably less than 2. 

According to our discussion in Section 5, for power law spectra s^2 in 

the thin-target model, and s^-3 in the thick-target model. For these spectral 
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parameters, from Figure 12 we see that the bulk of the positron emitters 
are rr + mesons, and from Figure 13 we get that the positron emitter-to- 
neutron ratio is about 0*2. The observed ratio, <j>g ^/cj^ 23^ then 

implies that f 7 f, i.e. the 0.51 MeV yield per positron is about the 
same as the 2.23 MeV yield per neutron. In view of the uncertainties 
involved in the deductions of f 7 , this is not an unreasonable result. 

We can also deduce f 7 for the various values of P obtained in Section 5, 
and in all cases we find acceptable values (0.1 < f 7 < 2) . For these 
exponential spectra the contribution of tt + mesons to the total positron 
production is greater than about 50% in all cases, except in the thick-target 
model with no He 3 in the photosphere where they contribute only about 25% 
of the positrons. A definite test for the possibility that the bulk of 
the positrons are due to tt + mesons would be the observation of gamma rays 
from T 7 ° decay. As there is no published data on high energy solar gamma 
rays, we shall not discuss this possibility in the present paper. 

As mentioned above, prompt 0.51 MeV photons could also result from 
positrons produced in e + -e pairs. Such pairs would be produced mainly 
by > 1 MeV electrons. Because these electrons also produce continuum 
emission by bremsstrahlung, it is possible to calculate the pair-to- 
bremsstrahlung yield (Bai and Ramaty 1975). By comparing these calculations 
with the observed continuum emission from the 1972, August 4 flare 
(Chupp 1974), we c ±nd that if the bremsstrahlung is produced by an isotropic 
distribution of electrons at the Sun, pair production could account for 
less than 1% of the positrons required to produce the 0.51 MeV line. 

However, if the electron beam is directed downward into the Sun, the 
bremsstrahlung efficiency in the backward direction is greatly reduced 
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f or relativistic particles (Petrosian 1973) . In this case pair production 
could account for essentially all the prompt 0.51 MeV photons of the 
1972, August 4 flare. 

Finally, we wish to mention that delayed photons (by about 15 minutes) 
were observed for the 1972, August 7 flare (Chupp et *1. 1973). These 
photons are most likely due to the long lived radioactive nuclei such 
C 11 and N 13 . 
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7. SUMMARY 

The observed gamma-ray lines from the 1972, August 4 flare, at 0.5, 

2.2, 4.4 and 6.1 MeV are due to positron annihilation, neutron capture 
on hydrogen, and deexcitation of excited states in C x2 a** , respec- 
tively. The strongest line is at 2.23 MeV. It is due to neutrons 

produced by nuclear reactions of flare accelerated particles with the 
ambient solar atmosphere. These neutrons are thermalized and captured 
by ambient protons in the photosphere to produce deuterons and 2.23 MeV 
gamma rays. Photospheric He 3 competes with the protons in capturing 
neutrons. Because captures on He 3 do not lead to photon emission, the 
observation of 2.23 MeV line emission from the Sun implies that the 
He 3 abundance in the photosphere cannot be much larger than that observed 
in the solar wind (He 3 /H ~ 5xl0~ 5 ). 

We have evaluated in detail the yield of neutrons and excited C 12 and 
0 16 nuclei from nuclear reactions of accelerated particles with the ambient, 
solar atmosphere. For the 1972, August 4 flare the neutrons are produced 
mainly in pa and ap reactions by primarv particles with energies greater 
than about 30 MeV/nucleon. The observed gamir * rays at 4.43 MeV and 6.14 MeV 
are principally due to proton induced interactions. Reactions induced by 
fast carbon and oxygen nuclei lead to Doppler broadened lines which cannot 
be distinguished from the continuum. 

From the comparison of the calculated and observed ratios of the lines 
at 4.43 MeV or 6.14 MeV to the line at 2.23 MeV it is possible to deduce the 
spectrum of the accelerated particles in the flare region. The spectral 
index s defined in equation (1) is 2 +0*2 for the thin-target model and 
3+0.3 for thick-target model if there is no He 3 in the photosphere; or 
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1.8 +0*2 and 2.7 +0.2 for these models if the ph^tospV ric lie 3 /H ratio 
is 5xl0~ 5 . The characteristic rigidity P Q defined in equation (2) is 
160 +35MV for the thin-t, et model and 110 +30MV f c r the thick-target 
model if He 3 /H = 0; or 230 +50MV and 180 +50MV for these models if 
He 3 /H = 5xl0' 5 . The total number of protons above 30 MeV relcrsed down- 
ard into the Sun in the thick-target model is about 10 33 . For the thin- 
target model about 2xl0 33 j ^otons escape from the fla^e region if the 

thickness of the target for relativistic paLticj.es is 1.5 g cm 2 as deduced 
from deuteron and helium -3 observations. 

The positrons which produce the 0.51 MeV line could be due to n + 
mesons, radioactive nuclei, and pair production. We have evaluated in 

detail the production of mesons and radioactive nuclei in nuclear reactions 

of accelerated charged particles with ambient nuclei. The lelativcly prompt 

nature of the 0.51 MeV line seems to favor positron production from rr + 
decay or e + -e pairs in the initial phase of the 1972, August 4 event. 

The delayed 0.51 MeV line emission observed for die 1972, . :gust 7 e^ent is very 

likely due to positrons from radioactive nuclei. 
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Figure Captions 


1. Neutron production cross sections 

2. Partial neutron production modes in the thin-target model with 

power -law spectra 

3. Partial neutron production modes in the thin- target model with 

exponential spectra 

4. Total neutron production in the thin and thick- target models with 

power-law and exponential spectra 

5. Probabilities for neutron escape, decay, and capture in the solar 

atmosphere (solid lines), and photon yields per neutron (dashed 
lines) for no He 3 in the photosphere 

6. Probabilities for neutron escape, decay, and capture in the solar 

atmosphere (solid line) , and photon yields per neutron (dashed 
lines) for He 3 /R * 5xl0~ 5 

7. C 12 * and 0 16 * production cross sections 

8. C 12 * and 0 16 * partial production modes in the thin-target model 

with power -law spectra 

9. Li 7 production cross section in at* reactions 

10. Ratios of the 4.43 MeV line intensity to the 2.23 MeV line intensity 

for thin and thick-target models, and power law and exponential 
spectra 

11. Partial positron production modes in the thin-target model with 

power-law spectra 

12. Ratios of the tt + meson yields to the total positron yields for the 

thin and thick- target models, and power law and exponential spectra 

13. Ratios of the total positron yield to the total neutron yield for 

the thin and thick-target models, and power law and exponential 
spectra 

14. Posicron free-annihilation rate and positronium formation rate in 

a hydrogen plasma 
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TABLE 2 


Prompt Gamma Ray Lines 


Photon Energy 
(MeV) 

Origin 

Production mode 

0.431 

7* 

Be deexc i ta t ion 

„ 4, , _ 7*0.4j1 

He (a,n) Be 

0.478 

7* 

Li deexci ta tion 

H. 4 (p,p> Li 7 * 0 ' 478 

1.63 

KT 20*1.63 . 

Ne deexcitation 

Ne^p.p'JNe 20 * 1,63 


14*3.94 -* M 14*2. 31, „14, , N „14*3.94 

N N deexc ltationN (p,p )N 

2.31 

14*2.31 

N deexc i ta tion 

N U (P. P ')» U * 2 - 31 

» 14 <p,p')n 14 * 3 - 94 V 4 * 2 

N 14 (p,n)0 14 - S 1412 - 31 

4.43 

12*4.43 - 

C deexci ta tion 

c 12 <p,pV 2 * 4 - 43 

-12, ,,12*4.43 

C (a,a')C 

16, , 12*4. 43 

0 (p, ) c 

5.2 

15*5.26 . 

0 deexci tat ion 

16, , 15*5.20 

0 (P, ) o 


15*5.28 - 

N deexcitation 

,16, ,15*5.28 

0 (p, )N 

6.14 

16*6 14 

0 * deexcitation 

-16, ,, 16*6.14 

0 (p,p') 0 

-16, ,, 16*6.14 

0 (a, a') 0 

7.12 

76*7.12 . _ _ 

0 deexcitati^n 

„16, ,,16*7.12 

0 (p,p')0 
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